This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.420354

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.12 while 0.88 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.289

0.271

0.099

0.578

1.019

0.126

0.143

NSEI

p_pelagic

13

1

2

0.116

0.100

0.032

0.287

1.027

0.068

0.234

NSEI

p_pelagic

13

10

1

0.304

0.289

0.113

0.581

1.011

0.121

0.124

NSEI

p_pelagic

13

10

2

0.123

0.108

0.034

0.287

1.015

0.067

0.224

NSEI

p_pelagic

13

11

2

0.128

0.113

0.036

0.300

1.014

0.070

0.214

NSEI

p_pelagic

13

12

1

0.312

0.298

0.113

0.588

1.010

0.126

0.111

NSEI

p_pelagic

13

12

2

0.129

0.114

0.038

0.306

1.012

0.071

0.212

NSEI

p_pelagic

13

13

1

0.316

0.301

0.115

0.580

1.015

0.124

0.121

NSEI

p_pelagic

13

14

2

0.130

0.117

0.039

0.292

1.012

0.066

0.198

NSEI

p_pelagic

13

16

1

0.317

0.303

0.125

0.581

1.016

0.118

0.118

NSEI

p_pelagic

13

17

1

0.325

0.313

0.136

0.581

1.013

0.116

0.103

NSEI

p_pelagic

13

18

1

0.338

0.329

0.140

0.595

1.010

0.117

0.077

NSEI

p_pelagic

13

19

2

0.149

0.137

0.052

0.314

1.012

0.069

0.173

NSEI

p_pelagic

13

2

1

0.292

0.276

0.101

0.579

1.014

0.124

0.129

NSEI

p_pelagic

13

2

2

0.117

0.102

0.033

0.290

1.011

0.066

0.226

NSEI

p_pelagic

13

21

2

0.169

0.156

0.057

0.353

1.010

0.076

0.170

NSEI

p_pelagic

13

3

2

0.118

0.102

0.032

0.298

1.026

0.069

0.233

NSEI

p_pelagic

13

4

2

0.120

0.104

0.033

0.300

1.016

0.070

0.227

NSEI

p_pelagic

13

5

1

0.292

0.272

0.104

0.585

1.012

0.125

0.160

NSEI

p_pelagic

13

5

2

0.119

0.103

0.033

0.284

1.019

0.067

0.240

NSEI

p_pelagic

13

6

1

0.294

0.280

0.103

0.576

1.015

0.126

0.111

NSEI

p_pelagic

13

6

2

0.120

0.104

0.032

0.298

1.019

0.068

0.234

NSEI

p_pelagic

13

7

1

0.295

0.277

0.103

0.577

1.014

0.126

0.143

NSEI

p_pelagic

13

7

2

0.121

0.106

0.033

0.297

1.014

0.068

0.220

NSEI

p_pelagic

13

8

2

0.121

0.106

0.034

0.289

1.022

0.069

0.218

NSEI

p_pelagic

13

9

1

0.305

0.291

0.107

0.576

1.013

0.125

0.112

NSEI

p_pelagic

13

9

2

0.123

0.109

0.036

0.296

1.013

0.067

0.208

NSEI

p_slope

14

13

2

0.073

0.054

0.009

0.255

1.013

0.064

0.298

NSEO

p_pelagic

16

1

1

0.353

0.342

0.115

0.640

1.059

0.138

0.080

SSEO

p_pelagic

16

1

2

0.330

0.314

0.112

0.614

1.060

0.136

0.117

SSEO

p_pelagic

16

10

1

0.358

0.349

0.134

0.637

1.049

0.132

0.068

SSEO

p_pelagic

16

10

2

0.337

0.323

0.120

0.627

1.039

0.133

0.105

SSEO

p_pelagic

16

11

1

0.358

0.347

0.132

0.641

1.049

0.135

0.082

SSEO

p_pelagic

16

11

2

0.342

0.329

0.122

0.630

1.046

0.138

0.095

SSEO

p_pelagic

16

12

1

0.363

0.356

0.135

0.641

1.046

0.134

0.052

SSEO

p_pelagic

16

12

2

0.337

0.326

0.118

0.613

1.033

0.131

0.084

SSEO

p_pelagic

16

13

1

0.360

0.348

0.135

0.658

1.056

0.137

0.087

SSEO

p_pelagic

16

13

2

0.339

0.327

0.121

0.622

1.048

0.132

0.091

SSEO

p_pelagic

16

14

1

0.352

0.343

0.132

0.617

1.033

0.131

0.069

SSEO

p_pelagic

16

14

2

0.345

0.332

0.131

0.619

1.042

0.131

0.099

SSEO

p_pelagic

16

15

1

0.355

0.344

0.143

0.631

1.029

0.128

0.086

SSEO

p_pelagic

16

15

2

0.347

0.333

0.129

0.616

1.031

0.131

0.107

SSEO

p_pelagic

16

16

1

0.353

0.341

0.139

0.629

1.028

0.129

0.093

SSEO

p_pelagic

16

16

2

0.351

0.336

0.139

0.623

1.025

0.129

0.116

SSEO

p_pelagic

16

17

1

0.415

0.406

0.188

0.677

1.033

0.128

0.070

SSEO

p_pelagic

16

17

2

0.367

0.358

0.145

0.637

1.034

0.128

0.070

SSEO

p_pelagic

16

18

1

0.392

0.385

0.174

0.657

1.028

0.129

0.054

SSEO

p_pelagic

16

18

2

0.368

0.359

0.150

0.633

1.033

0.127

0.071

SSEO

p_pelagic

16

19

1

0.382

0.372

0.169

0.633

1.021

0.124

0.080

SSEO

p_pelagic

16

19

2

0.366

0.359

0.150

0.631

1.024

0.125

0.056

SSEO

p_pelagic

16

2

1

0.352

0.343

0.116

0.643

1.057

0.142

0.063

SSEO

p_pelagic

16

2

2

0.330

0.314

0.110

0.613

1.043

0.134

0.119

SSEO

p_pelagic

16

20

1

0.388

0.380

0.176

0.649

1.016

0.123

0.065

SSEO

p_pelagic

16

20

2

0.368

0.358

0.155

0.636

1.019

0.126

0.080

SSEO

p_pelagic

16

21

2

0.370

0.362

0.159

0.623

1.018

0.123

0.065

SSEO

p_pelagic

16

3

1

0.352

0.344

0.112

0.642

1.055

0.142

0.056

SSEO

p_pelagic

16

3

2

0.333

0.321

0.111

0.625

1.048

0.135

0.089

SSEO

p_dsr

16

30

1

0.420

0.419

0.393

0.450

1.028

0.014

0.070

SSEO

p_dsr

16

31

1

0.402

0.401

0.374

0.433

1.013

0.016

0.064

SSEO

p_slope

16

36

1

0.061

0.062

0.041

0.079

1.036

0.009

-0.106

SSEO

p_slope

16

37

1

0.065

0.066

0.046

0.083

1.012

0.009

-0.110

SSEO

p_pelagic

16

4

1

0.353

0.339

0.122

0.656

1.052

0.143

0.098

SSEO

p_pelagic

16

4

2

0.335

0.323

0.111

0.636

1.047

0.137

0.087

SSEO

p_pelagic

16

42

1

0.833

0.833

0.826

0.840

1.055

0.003

0.000

SSEO

p_yellow

16

42

1

0.407

0.407

0.385

0.429

1.022

0.011

0.000

SSEO

p_dsr

16

42

1

0.445

0.438

0.406

0.503

2.015

0.027

0.260

SSEO

p_slope

16

42

1

0.174

0.183

0.115

0.209

2.222

0.028

-0.323

SSEO

p_pelagic

16

5

1

0.356

0.347

0.121

0.634

1.053

0.138

0.065

SSEO

p_pelagic

16

5

2

0.333

0.321

0.115

0.619

1.056

0.135

0.089

SSEO

p_pelagic

16

6

1

0.351

0.338

0.121

0.645

1.071

0.139

0.093

SSEO

p_pelagic

16

6

2

0.335

0.325

0.113

0.613

1.059

0.135

0.074

SSEO

p_pelagic

16

7

1

0.355

0.349

0.123

0.644

1.039

0.138

0.043

SSEO

p_pelagic

16

7

2

0.337

0.325

0.113

0.617

1.046

0.134

0.089

SSEO

p_pelagic

16

8

1

0.356

0.345

0.126

0.651

1.049

0.140

0.079

SSEO

p_pelagic

16

8

2

0.339

0.327

0.112

0.631

1.039

0.135

0.089

SSEO

p_pelagic

16

9

1

0.358

0.343

0.130

0.666

1.050

0.142

0.106

SSEO

p_pelagic

16

9

2

0.338

0.326

0.116

0.625

1.039

0.136

0.088

SSEO

p_pelagic

3

1

1

0.553

0.560

0.286

0.787

1.014

0.133

-0.053

PWSI

p_pelagic

3

1

2

0.187

0.175

0.061

0.380

1.012

0.084

0.142

PWSI

p_pelagic

3

10

1

0.549

0.555

0.275

0.784

1.013

0.133

-0.045

PWSI

p_pelagic

3

10

2

0.189

0.177

0.065

0.391

1.015

0.085

0.141

PWSI

p_pelagic

3

11

1

0.550

0.554

0.281

0.790

1.015

0.132

-0.030

PWSI

p_pelagic

3

12

1

0.551

0.560

0.282

0.785

1.016

0.131

-0.069

PWSI

p_pelagic

3

13

1

0.552

0.559

0.283

0.785

1.010

0.129

-0.054

PWSI

p_pelagic

3

16

2

0.188

0.175

0.067

0.377

1.018

0.082

0.159

PWSI

p_pelagic

3

17

1

0.547

0.550

0.289

0.784

1.012

0.128

-0.023

PWSI

p_pelagic

3

17

2

0.189

0.175

0.065

0.385

1.010

0.083

0.169

PWSI

p_pelagic

3

2

1

0.548

0.554

0.278

0.783

1.013

0.132

-0.045

PWSI

p_pelagic

3

3

1

0.548

0.554

0.271

0.790

1.017

0.135

-0.045

PWSI

p_pelagic

3

4

1

0.552

0.558

0.286

0.787

1.019

0.133

-0.045

PWSI

p_pelagic

3

4

2

0.191

0.178

0.061

0.404

1.013

0.087

0.149

PWSI

p_pelagic

3

5

1

0.553

0.560

0.287

0.781

1.012

0.131

-0.053

PWSI

p_pelagic

3

5

2

0.191

0.178

0.062

0.392

1.022

0.086

0.151

PWSI

p_pelagic

3

6

1

0.551

0.556

0.281

0.785

1.019

0.129

-0.039

PWSI

p_pelagic

3

6

2

0.189

0.175

0.065

0.395

1.014

0.085

0.165

PWSI

p_pelagic

3

7

1

0.549

0.554

0.282

0.781

1.011

0.132

-0.038

PWSI

p_pelagic

3

7

2

0.188

0.176

0.059

0.380

1.013

0.083

0.144

PWSI

p_pelagic

3

8

1

0.551

0.555

0.295

0.788

1.013

0.132

-0.030

PWSI

p_pelagic

3

8

2

0.188

0.176

0.062

0.391

1.010

0.085

0.141

PWSI

p_pelagic

3

9

1

0.553

0.558

0.276

0.796

1.018

0.134

-0.037

PWSI

p_pelagic

3

9

2

0.189

0.175

0.064

0.383

1.015

0.083

0.168

PWSI

p_pelagic

4

1

1

0.546

0.551

0.284

0.793

1.012

0.132

-0.038

PWSO

p_pelagic

4

10

1

0.549

0.552

0.292

0.789

1.013

0.130

-0.023

PWSO

p_pelagic

4

10

2

0.554

0.558

0.294

0.803

1.015

0.132

-0.030

PWSO

p_pelagic

4

11

1

0.550

0.557

0.290

0.780

1.015

0.130

-0.054

PWSO

p_pelagic

4

11

2

0.558

0.560

0.285

0.795

1.010

0.132

-0.015

PWSO

p_pelagic

4

12

1

0.546

0.549

0.288

0.781

1.015

0.129

-0.023

PWSO

p_pelagic

4

12

2

0.554

0.558

0.287

0.798

1.014

0.131

-0.031

PWSO

p_pelagic

4

13

1

0.545

0.551

0.283

0.789

1.013

0.131

-0.046

PWSO

p_pelagic

4

13

2

0.552

0.558

0.279

0.797

1.019

0.132

-0.045

PWSO

p_pelagic

4

15

2

0.551

0.556

0.293

0.784

1.015

0.129

-0.039

PWSO

p_pelagic

4

16

1

0.544

0.551

0.284

0.771

1.020

0.130

-0.054

PWSO

p_pelagic

4

16

2

0.556

0.559

0.303

0.793

1.017

0.132

-0.023

PWSO

p_pelagic

4

17

2

0.550

0.553

0.287

0.794

1.020

0.131

-0.023

PWSO

p_pelagic

4

18

1

0.546

0.552

0.296

0.789

1.014

0.130

-0.046

PWSO

p_pelagic

4

18

2

0.556

0.563

0.292

0.800

1.017

0.134

-0.052

PWSO

p_pelagic

4

19

2

0.551

0.551

0.295

0.791

1.012

0.131

0.000

PWSO

p_pelagic

4

2

1

0.549

0.552

0.286

0.785

1.014

0.130

-0.023

PWSO

p_pelagic

4

3

1

0.552

0.557

0.281

0.793

1.014

0.132

-0.038

PWSO

p_pelagic

4

3

2

0.552

0.556

0.288

0.796

1.020

0.133

-0.030

PWSO

p_pelagic

4

4

2

0.551

0.556

0.288

0.793

1.019

0.131

-0.038

PWSO

p_pelagic

4

5

1

0.543

0.547

0.285

0.783

1.016

0.132

-0.030

PWSO

p_pelagic

4

5

2

0.555

0.560

0.296

0.795

1.019

0.133

-0.038

PWSO

p_pelagic

4

6

1

0.552

0.560

0.290

0.790

1.010

0.133

-0.060

PWSO

p_pelagic

4

7

1

0.549

0.553

0.288

0.791

1.012

0.131

-0.031

PWSO

p_pelagic

4

7

2

0.555

0.559

0.293

0.790

1.018

0.131

-0.030

PWSO

p_pelagic

4

8

1

0.551

0.555

0.280

0.788

1.011

0.131

-0.031

PWSO

p_pelagic

4

8

2

0.552

0.557

0.280

0.791

1.013

0.135

-0.037

PWSO

p_pelagic

4

9

1

0.551

0.556

0.288

0.795

1.018

0.134

-0.037

PWSO

p_pelagic

4

9

2

0.552

0.556

0.290

0.797

1.013

0.133

-0.030

PWSO

p_black

8

40

1

0.475

0.415

0.341

0.963

1.016

0.178

0.337

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

1

1.012

0.670

0.287

4.242

1.193

1.047

0.327

northeast

beta0_pH

13

4

0.603

0.692

-0.964

1.676

1.141

0.664

-0.134

NSEI

beta1_pH

13

4

1.396

1.341

0.234

2.904

1.138

0.677

0.081

NSEI

beta0_pH

14

4

0.399

0.453

-0.762

1.286

1.129

0.530

-0.102

NSEO

beta1_pH

14

4

1.967

1.921

0.911

3.246

1.126

0.603

0.076

NSEO

beta0_pH

15

4

0.573

0.692

-1.026

1.548

1.314

0.650

-0.183

SSEI

beta1_pH

15

4

1.405

1.296

0.314

3.006

1.291

0.691

0.158

SSEI

beta0_yellow

15

-1.282

-1.237

-1.967

-1.027

1.254

0.259

-0.174

SSEI

beta1_yellow

15

0.192

0.073

0.011

1.081

1.156

0.297

0.401

SSEI

beta1_pH

16

4

3.529

3.520

2.216

4.964

1.107

0.699

0.013

SSEO

beta0_pelagic

16

-0.689

-0.627

-1.624

-0.013

1.161

0.454

-0.136

SSEO

beta1_pelagic

16

1.890

1.865

0.832

3.055

1.181

0.634

0.039

SSEO

mu_beta2_pH

2

1

0.946

0.588

0.061

3.962

1.115

1.006

0.356

NG

beta2_yellow

2

-5.744

-3.671

-20.442

-0.264

1.338

5.590

-0.371

NG

mu_beta0_pH

3

4

0.690

0.757

-0.543

1.622

1.136

0.559

-0.120

PWSI

tau_beta1_pH

3

4

5.778

2.866

0.231

34.183

1.137

9.370

0.311

PWSI

beta1_pelagic

3

0.769

0.718

0.309

1.587

1.108

0.301

0.169

PWSI

beta0_pelagic

4

0.219

0.231

-0.353

0.668

1.107

0.243

-0.049

PWSO

beta1_pelagic

4

0.832

0.822

0.334

1.506

1.112

0.268

0.037

PWSO

beta2_pH

8

1

0.683

0.344

0.123

4.197

1.132

1.045

0.324

afognak

beta2_pH

9

1

0.899

0.512

0.148

4.022

1.127

1.037

0.373

eastside

beta2_pH

9

3

0.282

0.195

0.088

0.865

1.125

0.520

0.167

eastside

lambda_H

15

0.019

0.014

0.003

0.061

1.226

0.043

0.117

SSEI

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

449.00

200.00

1,016.00

342.00

Prop. of total

0.69

0.31

0.75

0.25

Prop. median less than 1

0.27

0.26

0.35

0.45

Prop. median less than 10

0.42

0.56

0.58

0.76

Prop. median less than 50

0.60

0.81

0.75

0.85

Median of the median values

24.67

6.32

4.22

1.48

Mean of the median values

348.75

134.02

184.73

105.38

Mean skew of posteriors

0.26

0.28

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

12.00

14.00

62.00

47.00

Prop. of total

0.46

0.54

0.57

0.43

Prop. median less than 1

0.42

0.43

0.71

0.66

Prop. median less than 10

0.42

1.00

0.84

0.89

Prop. median less than 50

0.42

1.00

0.87

0.96

Median of the median values

180.53

1.02

0.06

0.33

Mean of the median values

421.52

2.29

37.09

4.61

Mean skew of posteriors

0.23

0.19

0.16

0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

1

beta0_pH

3

beta0_yellow

1

beta1_pelagic

2

beta1_pH

3

beta1_yellow

1

beta2_pH

4

beta2_yellow

1

eps_dsr

1

eps_pH

14

eps_sl

1

H_ay

4

H_ayg

1

H_ayu

4

Hb_ay

4

Hb_ayg

1

Hb_ayu

4

Hd_ayg

1

Ho_ayg

1

Ho_ayu

3

Hp_ay

4

Hp_ayg

1

Hp_ayu

4

Hs_ay

1

Hs_ayg

1

Hs_ayu

2

Hy_ay

4

Hy_ayg

1

Hy_ayu

8

lambda_H

1

mean_eps_pH

1

mu_beta0_pH

1

mu_beta2_pH

1

p_dsr

1

p_slope

1

pH

9

R_ay

4

R_ayg

2

R_ayu

5

Rb_ay

4

Rb_ayg

2

Rb_ayu

5

Rd_ay

2

Rd_ayg

3

Rd_ayu

3

re_dsr

1

re_pH

10

re_slope

1

Ro_ay

24

Ro_ayg

13

Ro_ayu

39

Rp_ay

4

Rp_ayg

2

Rp_ayu

4

Rs_ay

3

Rs_ayg

10

Rs_ayu

20

Ry_ay

10

Ry_ayg

3

Ry_ayu

25

tau_beta1_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.020 0.072 -0.149 -0.024 0.129
mu_bc_H[2] 0.075 0.053 -0.022 0.072 0.186
mu_bc_H[3] -0.266 0.068 -0.394 -0.268 -0.129
mu_bc_H[4] -0.931 0.183 -1.312 -0.926 -0.580
mu_bc_H[5] 1.355 1.197 0.129 1.082 4.330
mu_bc_H[6] -2.375 0.321 -2.994 -2.381 -1.750
mu_bc_H[7] -0.453 0.103 -0.664 -0.449 -0.261
mu_bc_H[8] 0.245 0.350 -0.334 0.207 1.023
mu_bc_H[9] -0.247 0.124 -0.487 -0.251 0.002
mu_bc_H[10] 0.060 0.065 -0.065 0.059 0.190
mu_bc_H[11] -0.128 0.037 -0.202 -0.128 -0.060
mu_bc_H[12] -0.256 0.111 -0.491 -0.251 -0.045
mu_bc_H[13] -0.162 0.072 -0.297 -0.163 -0.023
mu_bc_H[14] -0.306 0.092 -0.483 -0.304 -0.124
mu_bc_H[15] -0.336 0.046 -0.427 -0.337 -0.243
mu_bc_H[16] -0.176 0.366 -0.828 -0.202 0.574
mu_bc_R[1] 1.036 0.123 0.793 1.034 1.272
mu_bc_R[2] 1.011 0.075 0.859 1.013 1.149
mu_bc_R[3] 0.227 0.093 0.038 0.229 0.406
mu_bc_R[4] 0.662 0.199 0.261 0.670 1.024
mu_bc_R[5] 1.674 0.374 0.965 1.673 2.452
mu_bc_R[6] -1.412 0.338 -2.067 -1.416 -0.737
mu_bc_R[7] 0.511 0.170 0.171 0.515 0.837
mu_bc_R[8] 0.647 0.144 0.359 0.645 0.936
mu_bc_R[9] 0.385 0.171 0.023 0.394 0.698
mu_bc_R[10] 0.897 0.100 0.707 0.896 1.092
mu_bc_R[11] 1.110 0.074 0.967 1.109 1.262
mu_bc_R[12] 0.867 0.164 0.539 0.869 1.188
mu_bc_R[13] 1.032 0.091 0.849 1.033 1.208
mu_bc_R[14] 0.862 0.123 0.620 0.863 1.116
mu_bc_R[15] 0.804 0.097 0.614 0.804 0.992
mu_bc_R[16] 1.353 0.114 1.118 1.354 1.571
tau_pH[1] 2.453 0.175 2.124 2.452 2.806
tau_pH[2] 1.607 0.151 1.332 1.600 1.919
tau_pH[3] 1.193 0.108 0.992 1.189 1.411
tau_pH[4] 1.026 0.147 0.771 1.011 1.342
tau_pH[5] 0.422 0.075 0.295 0.414 0.583
beta0_pH[1,1] 0.551 0.166 0.212 0.556 0.868
beta0_pH[2,1] 1.428 0.164 1.063 1.438 1.727
beta0_pH[3,1] 1.410 0.208 0.928 1.434 1.765
beta0_pH[4,1] 1.595 0.219 1.125 1.605 1.998
beta0_pH[5,1] -0.819 0.407 -1.824 -0.777 -0.134
beta0_pH[6,1] -0.588 0.521 -2.107 -0.445 0.043
beta0_pH[7,1] 0.757 0.010 0.736 0.757 0.777
beta0_pH[8,1] -0.652 0.410 -1.776 -0.581 -0.092
beta0_pH[9,1] -0.341 0.312 -1.020 -0.327 0.208
beta0_pH[10,1] 0.245 0.205 -0.190 0.253 0.599
beta0_pH[11,1] 0.147 0.304 -0.548 0.171 0.662
beta0_pH[12,1] 0.627 0.186 0.250 0.632 0.973
beta0_pH[13,1] -0.657 0.144 -0.964 -0.648 -0.397
beta0_pH[14,1] 0.014 0.175 -0.339 0.015 0.352
beta0_pH[15,1] 0.072 0.175 -0.307 0.086 0.381
beta0_pH[16,1] -0.854 0.214 -1.265 -0.863 -0.435
beta0_pH[1,2] 2.564 0.085 2.407 2.564 2.737
beta0_pH[2,2] 3.079 0.022 3.038 3.079 3.121
beta0_pH[3,2] 3.626 0.041 3.548 3.625 3.706
beta0_pH[4,2] 2.936 0.019 2.899 2.936 2.975
beta0_pH[5,2] 4.569 1.529 2.745 4.252 8.360
beta0_pH[6,2] 3.324 0.157 3.020 3.318 3.633
beta0_pH[7,2] 2.032 0.051 1.934 2.031 2.131
beta0_pH[8,2] 2.824 0.065 2.700 2.823 2.954
beta0_pH[9,2] 3.954 0.210 3.580 3.945 4.397
beta0_pH[10,2] 3.997 0.127 3.753 3.994 4.260
beta0_pH[11,2] 1.885 0.009 1.867 1.885 1.903
beta0_pH[12,2] 1.792 0.068 1.664 1.788 1.927
beta0_pH[13,2] 2.496 0.026 2.445 2.496 2.547
beta0_pH[14,2] 1.803 0.024 1.756 1.803 1.850
beta0_pH[15,2] 2.545 0.017 2.513 2.545 2.578
beta0_pH[16,2] 2.766 0.025 2.718 2.766 2.816
beta0_pH[1,3] 1.071 0.395 0.180 1.143 1.661
beta0_pH[2,3] 2.126 0.017 2.091 2.127 2.160
beta0_pH[3,3] 2.692 0.024 2.644 2.692 2.739
beta0_pH[4,3] 2.981 0.035 2.913 2.980 3.046
beta0_pH[5,3] 3.502 2.411 0.897 2.921 10.112
beta0_pH[6,3] 1.363 0.376 0.731 1.328 2.183
beta0_pH[7,3] 0.859 0.050 0.760 0.858 0.956
beta0_pH[8,3] 0.179 0.611 -1.509 0.290 0.631
beta0_pH[9,3] -2.299 1.038 -4.636 -2.153 -0.577
beta0_pH[10,3] -0.497 0.834 -2.386 -0.411 0.861
beta0_pH[11,3] 1.088 0.550 -0.139 1.125 2.037
beta0_pH[12,3] 0.963 0.652 -0.466 0.975 2.169
beta0_pH[13,3] 0.983 0.533 -0.249 1.021 1.941
beta0_pH[14,3] 0.563 0.464 -0.482 0.609 1.357
beta0_pH[15,3] 0.722 0.443 -0.240 0.741 1.500
beta0_pH[16,3] 1.178 0.634 -0.053 1.170 2.440
beta0_pH[11,4] 0.859 0.548 -0.408 0.908 1.826
beta0_pH[12,4] 0.905 0.692 -0.503 0.923 2.220
beta0_pH[13,4] 0.603 0.664 -0.964 0.692 1.676
beta0_pH[14,4] 0.399 0.530 -0.762 0.453 1.286
beta0_pH[15,4] 0.573 0.650 -1.026 0.692 1.548
beta0_pH[16,4] 0.862 0.621 -0.477 0.892 2.026
beta0_pH[11,5] 2.964 0.145 2.685 2.964 3.258
beta0_pH[12,5] -0.919 0.509 -1.966 -0.902 0.025
beta0_pH[13,5] 4.922 0.426 4.165 4.906 5.832
beta0_pH[14,5] 3.911 0.374 3.226 3.888 4.703
beta0_pH[15,5] 2.289 0.129 2.052 2.286 2.554
beta0_pH[16,5] 0.244 0.063 0.158 0.232 0.401
beta1_pH[1,1] 2.573 0.448 1.819 2.533 3.589
beta1_pH[2,1] 1.577 0.348 1.035 1.537 2.535
beta1_pH[3,1] 1.782 0.391 1.135 1.743 2.704
beta1_pH[4,1] 1.870 0.407 1.197 1.826 2.863
beta1_pH[5,1] 2.359 0.523 1.517 2.299 3.671
beta1_pH[6,1] 2.270 0.875 1.111 2.039 4.440
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.451 0.871 1.217 2.240 4.469
beta1_pH[9,1] 2.121 0.698 1.160 1.985 4.245
beta1_pH[10,1] 2.213 0.372 1.575 2.190 3.035
beta1_pH[11,1] 3.266 0.518 2.447 3.179 4.489
beta1_pH[12,1] 2.900 0.274 2.369 2.902 3.454
beta1_pH[13,1] 3.457 0.237 3.014 3.448 3.955
beta1_pH[14,1] 2.834 0.321 2.227 2.827 3.496
beta1_pH[15,1] 2.383 0.260 1.901 2.369 2.920
beta1_pH[16,1] 4.004 0.297 3.417 4.023 4.562
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.950 1.246 0.064 1.780 5.349
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.773 2.577 0.790 4.279 11.532
beta1_pH[9,3] 7.651 2.637 3.637 7.243 13.594
beta1_pH[10,3] 6.589 2.670 2.732 6.127 12.618
beta1_pH[11,3] 2.201 0.597 1.120 2.165 3.493
beta1_pH[12,3] 3.250 0.674 2.031 3.225 4.741
beta1_pH[13,3] 1.992 0.566 0.920 1.971 3.293
beta1_pH[14,3] 2.210 0.518 1.298 2.178 3.364
beta1_pH[15,3] 2.102 0.513 1.164 2.082 3.190
beta1_pH[16,3] 3.705 0.689 2.291 3.715 5.009
beta1_pH[11,4] 1.923 0.593 0.820 1.904 3.193
beta1_pH[12,4] 3.178 0.711 1.805 3.157 4.618
beta1_pH[13,4] 1.396 0.677 0.234 1.341 2.904
beta1_pH[14,4] 1.967 0.603 0.911 1.921 3.246
beta1_pH[15,4] 1.405 0.691 0.314 1.296 3.006
beta1_pH[16,4] 3.529 0.699 2.216 3.520 4.964
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.255 0.851 0.360 0.931 3.384
beta2_pH[2,1] 1.631 0.968 0.239 1.473 3.427
beta2_pH[3,1] 1.412 0.837 0.229 1.270 2.941
beta2_pH[4,1] 1.184 0.790 0.229 0.941 2.917
beta2_pH[5,1] 1.879 1.612 0.260 1.425 6.419
beta2_pH[6,1] 0.879 1.375 0.106 0.367 5.103
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.683 1.045 0.123 0.344 4.197
beta2_pH[9,1] 0.899 1.037 0.148 0.512 4.022
beta2_pH[10,1] 1.012 1.047 0.287 0.670 4.242
beta2_pH[11,1] 0.463 0.264 0.188 0.400 1.104
beta2_pH[12,1] 1.556 0.792 0.647 1.301 3.632
beta2_pH[13,1] 1.002 0.388 0.550 0.923 2.065
beta2_pH[14,1] 1.144 0.671 0.435 0.945 3.234
beta2_pH[15,1] 1.052 0.632 0.374 0.877 3.024
beta2_pH[16,1] 1.052 0.310 0.644 0.995 1.773
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 6.870 6.654 0.173 5.144 22.973
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 1.449 1.367 0.070 1.020 5.630
beta2_pH[9,3] 0.282 0.520 0.088 0.195 0.865
beta2_pH[10,3] 0.406 0.604 0.105 0.260 1.699
beta2_pH[11,3] -6.178 4.915 -19.394 -4.832 -0.583
beta2_pH[12,3] -6.983 4.798 -19.422 -5.802 -1.229
beta2_pH[13,3] -6.297 4.833 -19.064 -4.843 -0.959
beta2_pH[14,3] -6.487 4.860 -20.336 -5.024 -1.217
beta2_pH[15,3] -6.575 4.828 -19.335 -5.232 -1.100
beta2_pH[16,3] -6.495 4.745 -19.410 -5.048 -1.479
beta2_pH[11,4] -5.876 4.424 -17.149 -4.717 -0.776
beta2_pH[12,4] -6.192 4.367 -16.798 -5.033 -1.059
beta2_pH[13,4] -6.010 4.393 -16.858 -4.762 -0.901
beta2_pH[14,4] -6.081 4.249 -16.784 -4.973 -1.074
beta2_pH[15,4] -5.845 4.326 -15.949 -4.741 -0.802
beta2_pH[16,4] -6.169 4.126 -16.574 -4.996 -1.396
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.940 1.110 32.887 34.923 37.356
beta3_pH[2,1] 32.630 1.471 30.094 32.514 36.108
beta3_pH[3,1] 33.069 1.286 30.524 33.022 35.758
beta3_pH[4,1] 32.730 1.763 29.737 32.639 36.619
beta3_pH[5,1] 27.643 1.260 25.948 27.473 31.207
beta3_pH[6,1] 33.229 3.321 26.351 33.062 40.720
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.151 3.803 27.261 35.274 42.105
beta3_pH[9,1] 32.403 2.827 28.073 31.939 39.765
beta3_pH[10,1] 32.636 1.271 30.330 32.587 35.336
beta3_pH[11,1] 31.566 1.110 29.453 31.566 33.944
beta3_pH[12,1] 30.787 0.483 29.784 30.810 31.728
beta3_pH[13,1] 31.415 0.548 30.361 31.398 32.489
beta3_pH[14,1] 32.483 0.782 31.188 32.411 34.148
beta3_pH[15,1] 31.013 0.727 29.521 31.015 32.439
beta3_pH[16,1] 30.145 0.345 29.443 30.141 30.849
beta3_pH[1,2] 34.608 5.498 25.405 34.809 43.471
beta3_pH[2,2] 34.507 5.560 25.468 34.469 43.540
beta3_pH[3,2] 30.573 7.150 18.644 30.771 42.310
beta3_pH[4,2] 30.524 7.141 18.755 30.697 42.315
beta3_pH[5,2] 30.499 7.185 18.622 30.687 42.430
beta3_pH[6,2] 30.462 7.226 18.648 30.435 42.348
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.072 3.398 29.390 39.357 42.795
beta3_pH[2,3] 34.487 5.421 25.412 34.473 43.573
beta3_pH[3,3] 30.429 7.129 18.738 30.273 42.360
beta3_pH[4,3] 30.641 7.205 18.600 30.820 42.284
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 41.959 1.467 37.575 42.056 43.886
beta3_pH[9,3] 36.184 1.900 34.525 35.437 41.476
beta3_pH[10,3] 37.979 2.548 34.668 37.339 43.474
beta3_pH[11,3] 40.750 1.050 38.895 40.772 42.705
beta3_pH[12,3] 42.039 0.712 40.322 42.285 42.832
beta3_pH[13,3] 41.470 0.752 39.969 41.481 42.876
beta3_pH[14,3] 40.515 0.697 39.083 40.554 41.786
beta3_pH[15,3] 40.646 0.642 39.199 40.670 41.886
beta3_pH[16,3] 41.840 0.392 41.034 41.897 42.608
beta3_pH[11,4] 40.655 1.112 38.478 40.581 42.994
beta3_pH[12,4] 41.806 0.819 40.228 42.080 42.932
beta3_pH[13,4] 41.881 1.233 39.037 41.984 43.807
beta3_pH[14,4] 40.440 0.767 38.840 40.470 41.839
beta3_pH[15,4] 41.405 1.204 38.620 41.398 43.605
beta3_pH[16,4] 41.655 0.428 40.906 41.663 42.559
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.264 0.020 0.224 0.264 0.302
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.159 0.173 0.187
beta4_pH[4,1] 0.293 0.011 0.272 0.294 0.315
beta4_pH[5,1] 0.426 0.266 0.091 0.386 0.948
beta4_pH[6,1] 0.779 0.368 0.310 0.729 1.480
beta4_pH[7,1] 0.855 0.348 0.331 0.808 1.671
beta4_pH[8,1] 0.157 0.069 0.044 0.150 0.315
beta4_pH[9,1] 0.213 0.086 0.060 0.209 0.390
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.447 0.012 0.424 0.447 0.469
beta4_pH[12,1] 0.642 0.022 0.599 0.641 0.686
beta4_pH[13,1] 0.631 0.013 0.605 0.630 0.657
beta4_pH[14,1] 0.670 0.032 0.610 0.669 0.735
beta4_pH[15,1] 0.798 0.021 0.758 0.797 0.838
beta4_pH[16,1] 1.220 0.043 1.141 1.220 1.307
beta4_pH[1,2] 0.309 0.096 0.137 0.304 0.511
beta4_pH[2,2] 0.515 0.023 0.472 0.514 0.561
beta4_pH[3,2] 0.446 0.012 0.423 0.445 0.468
beta4_pH[4,2] 0.654 0.038 0.582 0.654 0.730
beta4_pH[5,2] 0.738 0.953 0.112 0.535 2.746
beta4_pH[6,2] 0.928 0.828 0.210 0.745 2.762
beta4_pH[7,2] 0.900 0.445 0.335 0.818 1.968
beta4_pH[8,2] 0.843 0.201 0.480 0.833 1.269
beta4_pH[9,2] 1.161 1.053 0.302 0.891 3.940
beta4_pH[10,2] 0.455 0.052 0.356 0.454 0.560
beta4_pH[11,2] 1.267 0.050 1.171 1.268 1.364
beta4_pH[12,2] 0.584 0.110 0.390 0.579 0.806
beta4_pH[13,2] 0.945 0.044 0.857 0.945 1.033
beta4_pH[14,2] 0.931 0.062 0.817 0.930 1.054
beta4_pH[15,2] 1.011 0.026 0.961 1.011 1.062
beta4_pH[16,2] 1.406 0.075 1.265 1.405 1.551
beta4_pH[1,3] 0.284 0.086 0.131 0.281 0.464
beta4_pH[2,3] 0.051 0.017 0.019 0.050 0.088
beta4_pH[3,3] 0.052 0.010 0.034 0.052 0.071
beta4_pH[4,3] 0.173 0.032 0.113 0.173 0.235
beta4_pH[5,3] 0.747 0.774 0.123 0.570 2.473
beta4_pH[6,3] 0.885 0.827 0.212 0.717 2.549
beta4_pH[7,3] 0.900 0.463 0.343 0.811 1.956
beta4_pH[8,3] 0.415 0.200 0.110 0.388 0.867
beta4_pH[9,3] 0.930 0.908 0.158 0.708 3.078
beta4_pH[10,3] 0.061 0.042 0.004 0.056 0.158
beta4_pH[11,3] 0.364 0.016 0.333 0.364 0.397
beta4_pH[12,3] 0.440 0.041 0.363 0.438 0.525
beta4_pH[13,3] 0.840 0.037 0.767 0.840 0.912
beta4_pH[14,3] 0.825 0.068 0.700 0.825 0.962
beta4_pH[15,3] 0.635 0.017 0.603 0.635 0.668
beta4_pH[16,3] 0.638 0.036 0.569 0.638 0.710
beta4_pH[1,4] 0.630 0.336 0.220 0.559 1.442
beta4_pH[2,4] 0.803 0.958 0.119 0.580 2.902
beta4_pH[3,4] 0.835 0.861 0.109 0.587 2.993
beta4_pH[4,4] 0.730 0.487 0.218 0.624 1.909
beta4_pH[5,4] 0.775 0.924 0.122 0.578 2.791
beta4_pH[6,4] 0.902 0.785 0.210 0.728 2.618
beta4_pH[7,4] 0.897 0.409 0.350 0.827 1.879
beta4_pH[8,4] 0.753 0.611 0.154 0.600 2.232
beta4_pH[9,4] 0.770 0.807 0.129 0.567 2.650
beta4_pH[10,4] 1.470 3.035 0.008 0.394 10.910
beta4_pH[11,4] 0.444 0.017 0.411 0.444 0.478
beta4_pH[12,4] 0.561 0.055 0.457 0.560 0.674
beta4_pH[13,4] 0.925 0.036 0.855 0.926 0.995
beta4_pH[14,4] 0.922 0.074 0.778 0.922 1.071
beta4_pH[15,4] 0.714 0.020 0.675 0.714 0.755
beta4_pH[16,4] 0.689 0.042 0.610 0.689 0.776
beta4_pH[1,5] 0.628 0.343 0.209 0.566 1.419
beta4_pH[2,5] 0.795 0.801 0.116 0.574 2.852
beta4_pH[3,5] 0.856 0.931 0.112 0.596 3.187
beta4_pH[4,5] 0.727 0.529 0.213 0.622 1.797
beta4_pH[5,5] 0.770 0.826 0.114 0.561 2.755
beta4_pH[6,5] 0.907 0.826 0.210 0.723 2.758
beta4_pH[7,5] 0.889 0.433 0.348 0.814 1.895
beta4_pH[8,5] 0.763 0.864 0.157 0.589 2.402
beta4_pH[9,5] 0.790 0.984 0.120 0.573 2.699
beta4_pH[10,5] 1.488 3.110 0.007 0.402 10.875
beta4_pH[11,5] 0.602 0.052 0.502 0.601 0.707
beta4_pH[12,5] 0.580 0.164 0.308 0.567 0.955
beta4_pH[13,5] 0.773 0.063 0.655 0.771 0.907
beta4_pH[14,5] 0.846 0.119 0.635 0.838 1.106
beta4_pH[15,5] 1.144 0.076 1.005 1.141 1.295
beta4_pH[16,5] 19.106 3.875 11.029 19.646 24.729
beta0_pelagic[1] 2.267 0.036 2.200 2.266 2.338
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.552
beta0_pelagic[3] 0.186 0.259 -0.446 0.244 0.563
beta0_pelagic[4] 0.219 0.243 -0.353 0.231 0.668
beta0_pelagic[5] 1.448 0.154 1.158 1.445 1.764
beta0_pelagic[6] 1.353 0.060 1.240 1.351 1.475
beta0_pelagic[7] 1.480 0.039 1.407 1.479 1.557
beta0_pelagic[8] 1.673 0.044 1.593 1.673 1.762
beta0_pelagic[9] 2.794 0.066 2.667 2.793 2.928
beta0_pelagic[10] 2.603 0.031 2.548 2.601 2.666
beta0_pelagic[11] -0.098 0.141 -0.400 -0.092 0.164
beta0_pelagic[12] 1.801 0.035 1.738 1.799 1.873
beta0_pelagic[13] -1.061 0.416 -1.711 -1.096 -0.202
beta0_pelagic[14] -0.280 0.168 -0.644 -0.269 0.018
beta0_pelagic[15] -0.211 0.034 -0.298 -0.197 -0.180
beta0_pelagic[16] -0.689 0.454 -1.624 -0.627 -0.013
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.769 0.301 0.309 0.718 1.587
beta1_pelagic[4] 0.832 0.268 0.334 0.822 1.506
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.104 0.220 0.731 1.082 1.641
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 2.885 0.608 1.648 2.860 3.915
beta1_pelagic[14] 1.758 0.434 1.091 1.671 2.823
beta1_pelagic[15] 0.094 0.161 0.000 0.016 0.609
beta1_pelagic[16] 1.890 0.634 0.832 1.865 3.055
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.139 5.921 0.169 4.433 20.583
beta2_pelagic[4] 7.403 6.177 0.317 5.985 22.792
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.881 1.526 0.284 1.389 5.920
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.175 0.070 0.095 0.159 0.358
beta2_pelagic[14] 0.924 0.950 0.217 0.616 3.624
beta2_pelagic[15] 2.665 1.901 0.165 2.308 7.240
beta2_pelagic[16] 0.534 0.793 0.112 0.261 3.174
beta3_pelagic[1] 28.851 3.860 22.298 28.939 35.120
beta3_pelagic[2] 28.888 3.869 22.290 28.935 35.145
beta3_pelagic[3] 27.807 2.664 23.202 28.094 32.028
beta3_pelagic[4] 26.163 1.529 23.060 26.280 29.848
beta3_pelagic[5] 29.449 6.552 18.608 29.325 40.396
beta3_pelagic[6] 29.476 6.606 18.627 29.473 40.316
beta3_pelagic[7] 29.394 6.617 18.474 29.332 40.333
beta3_pelagic[8] 29.577 6.553 18.602 29.716 40.402
beta3_pelagic[9] 29.481 6.624 18.627 29.562 40.443
beta3_pelagic[10] 29.421 6.700 18.496 29.442 40.391
beta3_pelagic[11] 30.314 1.332 28.147 30.122 33.640
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 30.974 3.763 24.769 30.670 38.563
beta3_pelagic[14] 34.085 1.599 31.470 33.939 37.889
beta3_pelagic[15] 34.127 5.770 23.478 34.954 41.766
beta3_pelagic[16] 30.459 3.023 24.924 30.533 37.200
beta0_yellow[1] -0.483 0.078 -0.636 -0.484 -0.330
beta0_yellow[2] 0.527 0.099 0.284 0.538 0.677
beta0_yellow[3] -0.218 0.137 -0.502 -0.218 0.012
beta0_yellow[4] 0.920 0.219 0.387 0.982 1.192
beta0_yellow[5] 0.256 0.328 -0.412 0.271 0.845
beta0_yellow[6] 0.818 0.107 0.617 0.818 1.028
beta0_yellow[7] 0.994 0.091 0.811 0.994 1.170
beta0_yellow[8] 1.025 0.092 0.848 1.022 1.209
beta0_yellow[9] 0.460 0.102 0.262 0.461 0.663
beta0_yellow[10] 0.539 0.071 0.409 0.537 0.679
beta0_yellow[11] -0.209 0.132 -0.488 -0.197 0.005
beta0_yellow[12] -1.081 0.095 -1.269 -1.081 -0.902
beta0_yellow[13] -1.339 0.174 -1.736 -1.323 -1.048
beta0_yellow[14] -0.252 0.151 -0.588 -0.238 -0.005
beta0_yellow[15] -1.282 0.259 -1.967 -1.237 -1.027
beta0_yellow[16] -0.250 0.093 -0.432 -0.250 -0.069
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.038 0.255 0.600 1.022 1.612
beta1_yellow[3] 0.488 0.268 0.021 0.494 0.964
beta1_yellow[4] 1.297 0.545 0.631 1.171 3.066
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.332 0.386 0.001 0.172 1.291
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.797 1.631 0.034 1.065 5.837
beta1_yellow[14] 0.703 0.960 0.012 0.431 4.044
beta1_yellow[15] 0.192 0.297 0.011 0.073 1.081
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -5.744 5.590 -20.442 -3.671 -0.264
beta2_yellow[3] -7.520 5.782 -20.913 -6.174 -0.333
beta2_yellow[4] -5.843 5.956 -20.771 -3.991 -0.135
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.772 5.828 -21.717 -6.512 -0.433
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.446 5.200 -18.678 -3.729 -0.289
beta2_yellow[14] -6.933 5.407 -19.533 -5.706 -0.334
beta2_yellow[15] -7.896 5.875 -21.184 -6.693 -0.314
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.985 1.330 27.128 28.750 32.455
beta3_yellow[3] 32.338 3.366 22.635 32.552 38.550
beta3_yellow[4] 29.668 3.177 23.776 28.500 35.538
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 34.749 3.556 30.356 33.768 43.303
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.493 2.665 29.113 32.383 38.039
beta3_yellow[14] 33.922 3.518 29.196 33.528 42.868
beta3_yellow[15] 36.583 3.979 30.376 36.289 43.683
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.014 0.076 -0.130 0.013 0.167
beta0_black[2] 1.855 0.032 1.792 1.854 1.917
beta0_black[3] 1.239 0.037 1.166 1.238 1.311
beta0_black[4] 2.425 0.048 2.334 2.424 2.520
beta0_black[5] 5.264 2.441 1.910 4.806 11.509
beta0_black[6] 5.314 2.188 2.708 4.804 11.104
beta0_black[7] 4.245 2.266 1.684 3.589 10.425
beta0_black[8] 0.879 0.211 0.488 0.870 1.315
beta0_black[9] 2.818 0.276 2.295 2.812 3.367
beta0_black[10] 1.287 0.051 1.188 1.287 1.388
beta0_black[11] 2.978 0.052 2.877 2.977 3.080
beta0_black[12] 4.914 0.193 4.536 4.913 5.293
beta0_black[13] 0.417 0.318 0.014 0.346 1.200
beta0_black[14] 2.853 0.105 2.649 2.851 3.056
beta0_black[15] 0.723 0.052 0.621 0.723 0.826
beta0_black[16] 3.772 0.074 3.624 3.772 3.910
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.573 1.841 -7.189 -2.118 -0.330
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.571 2.982 29.448 38.193 40.853
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.328 0.094 -0.511 -0.326 -0.145
beta4_black[2] 0.256 0.054 0.150 0.256 0.359
beta4_black[3] -0.719 0.064 -0.844 -0.719 -0.594
beta4_black[4] 0.629 0.128 0.388 0.630 0.885
beta4_black[5] 0.781 1.663 -1.173 0.491 4.373
beta4_black[6] 0.860 1.555 -0.883 0.532 4.653
beta4_black[7] 0.747 1.258 -0.898 0.509 3.912
beta4_black[8] 0.088 0.281 -0.508 0.102 0.606
beta4_black[9] 1.039 0.821 -0.082 0.854 3.067
beta4_black[10] 0.074 0.067 -0.059 0.074 0.204
beta4_black[11] -0.608 0.065 -0.741 -0.607 -0.485
beta4_black[12] 1.237 0.294 0.689 1.232 1.833
beta4_black[13] -1.001 0.076 -1.149 -1.000 -0.852
beta4_black[14] 0.224 0.176 -0.113 0.223 0.573
beta4_black[15] 0.002 0.047 -0.090 0.002 0.096
beta4_black[16] -0.415 0.118 -0.644 -0.414 -0.181
beta0_dsr[11] -0.329 0.009 -0.347 -0.329 -0.311
beta0_dsr[12] 1.409 0.045 1.324 1.408 1.496
beta0_dsr[13] 0.155 0.021 0.114 0.154 0.197
beta0_dsr[14] -0.127 0.017 -0.159 -0.127 -0.093
beta0_dsr[15] 0.206 0.014 0.179 0.206 0.233
beta0_dsr[16] 0.159 0.016 0.128 0.159 0.189
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.544 0.035 0.472 0.543 0.613
beta4_dsr[12] 0.799 0.227 0.379 0.793 1.246
beta4_dsr[13] -0.077 0.041 -0.157 -0.078 0.003
beta4_dsr[14] -0.208 0.070 -0.345 -0.208 -0.074
beta4_dsr[15] 0.217 0.022 0.174 0.217 0.258
beta4_dsr[16] -0.049 0.046 -0.138 -0.051 0.040
beta0_slope[11] -2.209 0.021 -2.249 -2.209 -2.167
beta0_slope[12] -4.985 0.184 -5.368 -4.974 -4.642
beta0_slope[13] -1.501 0.034 -1.571 -1.500 -1.435
beta0_slope[14] -2.712 0.051 -2.814 -2.710 -2.613
beta0_slope[15] -1.592 0.022 -1.637 -1.592 -1.550
beta0_slope[16] -2.647 0.041 -2.726 -2.648 -2.564
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.333 0.076 -0.487 -0.330 -0.186
beta4_slope[12] -4.667 0.805 -6.337 -4.609 -3.198
beta4_slope[13] 0.326 0.054 0.223 0.326 0.432
beta4_slope[14] -0.161 0.152 -0.475 -0.156 0.125
beta4_slope[15] -0.359 0.033 -0.425 -0.359 -0.294
beta4_slope[16] -0.147 0.102 -0.350 -0.148 0.045
beta5_slope[11] 3.068 0.109 2.866 3.063 3.294
beta5_slope[12] 2.002 0.388 1.287 1.985 2.788
beta5_slope[13] 0.885 0.078 0.731 0.886 1.040
beta5_slope[14] 3.749 0.358 3.109 3.730 4.498
beta5_slope[15] 1.307 0.072 1.167 1.307 1.447
beta5_slope[16] 2.730 0.233 2.290 2.728 3.192
tau_comp[1] 3.770 0.297 3.206 3.760 4.369
tau_comp[2] 2.853 0.250 2.383 2.842 3.387
tau_comp[3] 0.868 0.074 0.727 0.865 1.014
tau_comp[4] 1.786 0.274 1.303 1.765 2.372
tau_comp[5] 1.273 0.186 0.943 1.262 1.676
sigma_H[1] 0.206 0.049 0.121 0.202 0.311
sigma_H[2] 0.176 0.029 0.128 0.174 0.244
sigma_H[3] 0.189 0.041 0.118 0.186 0.278
sigma_H[4] 0.461 0.081 0.323 0.452 0.648
sigma_H[5] 0.959 0.201 0.589 0.948 1.405
sigma_H[6] 0.526 0.194 0.159 0.519 0.931
sigma_H[7] 0.310 0.064 0.211 0.299 0.458
sigma_H[8] 0.418 0.094 0.255 0.409 0.625
sigma_H[9] 0.514 0.119 0.329 0.495 0.791
sigma_H[10] 0.210 0.039 0.142 0.207 0.298
sigma_H[11] 0.267 0.045 0.192 0.262 0.366
sigma_H[12] 0.379 0.147 0.196 0.331 0.712
sigma_H[13] 0.219 0.037 0.158 0.215 0.301
sigma_H[14] 0.495 0.089 0.339 0.489 0.689
sigma_H[15] 0.232 0.038 0.165 0.228 0.317
sigma_H[16] 0.217 0.041 0.147 0.213 0.311
lambda_H[1] 2.301 3.508 0.088 1.056 11.787
lambda_H[2] 5.219 5.195 0.513 3.833 17.868
lambda_H[3] 5.630 8.004 0.254 2.899 28.466
lambda_H[4] 0.007 0.005 0.001 0.005 0.018
lambda_H[5] 2.474 5.454 0.029 0.726 17.861
lambda_H[6] 5.695 12.543 0.008 0.206 43.119
lambda_H[7] 0.012 0.009 0.002 0.010 0.034
lambda_H[8] 7.478 9.456 0.101 4.262 33.928
lambda_H[9] 0.015 0.010 0.003 0.013 0.041
lambda_H[10] 0.226 0.235 0.027 0.156 0.788
lambda_H[11] 0.246 0.342 0.011 0.117 1.188
lambda_H[12] 3.492 4.367 0.161 2.144 15.030
lambda_H[13] 3.225 2.846 0.300 2.426 10.891
lambda_H[14] 3.438 4.209 0.213 2.192 14.623
lambda_H[15] 0.019 0.043 0.003 0.014 0.061
lambda_H[16] 0.610 0.917 0.035 0.301 3.206
mu_lambda_H[1] 4.103 1.875 1.074 3.869 8.167
mu_lambda_H[2] 3.617 1.968 0.478 3.459 7.846
mu_lambda_H[3] 3.386 1.822 0.731 3.091 7.531
sigma_lambda_H[1] 8.137 4.309 1.877 7.346 18.085
sigma_lambda_H[2] 7.987 4.770 0.822 7.319 18.497
sigma_lambda_H[3] 6.150 3.990 0.941 5.235 16.154
beta_H[1,1] 6.543 1.270 3.591 6.767 8.411
beta_H[2,1] 9.758 0.613 8.371 9.796 10.852
beta_H[3,1] 7.822 0.778 6.016 7.924 9.097
beta_H[4,1] 10.472 7.579 -4.911 10.522 24.769
beta_H[5,1] -0.254 2.531 -5.618 -0.116 4.070
beta_H[6,1] 3.007 4.445 -7.917 4.505 8.473
beta_H[7,1] 1.135 5.683 -10.912 1.548 11.700
beta_H[8,1] 1.257 3.358 -2.577 1.255 3.439
beta_H[9,1] 13.310 5.748 2.282 13.190 25.485
beta_H[10,1] 6.763 1.824 2.843 6.875 10.199
beta_H[11,1] 5.065 3.435 -2.670 5.736 9.909
beta_H[12,1] 2.615 1.100 0.620 2.563 5.012
beta_H[13,1] 9.126 0.862 7.233 9.178 10.612
beta_H[14,1] 2.143 1.045 0.135 2.146 4.268
beta_H[15,1] -7.040 3.432 -13.241 -7.247 0.171
beta_H[16,1] 3.782 2.890 -0.662 3.375 10.815
beta_H[1,2] 7.784 0.262 7.254 7.791 8.286
beta_H[2,2] 9.882 0.164 9.572 9.881 10.199
beta_H[3,2] 8.796 0.203 8.389 8.798 9.201
beta_H[4,2] 3.068 1.428 0.413 3.043 5.999
beta_H[5,2] 1.778 0.988 -0.215 1.773 3.749
beta_H[6,2] 5.995 1.110 3.550 6.148 7.774
beta_H[7,2] 2.397 1.053 0.457 2.353 4.585
beta_H[8,2] 3.037 1.008 1.385 3.162 4.260
beta_H[9,2] 3.243 1.099 1.082 3.201 5.548
beta_H[10,2] 8.046 0.354 7.330 8.063 8.745
beta_H[11,2] 9.793 0.622 8.840 9.683 11.169
beta_H[12,2] 3.951 0.376 3.261 3.931 4.750
beta_H[13,2] 9.156 0.238 8.713 9.148 9.667
beta_H[14,2] 4.021 0.346 3.364 4.026 4.689
beta_H[15,2] 11.557 0.616 10.275 11.584 12.683
beta_H[16,2] 4.481 0.819 2.998 4.463 6.141
beta_H[1,3] 8.449 0.264 7.987 8.439 8.985
beta_H[2,3] 9.859 0.149 9.565 9.857 10.170
beta_H[3,3] 9.455 0.162 9.126 9.453 9.782
beta_H[4,3] -2.174 0.865 -3.830 -2.183 -0.429
beta_H[5,3] 3.882 0.595 2.673 3.912 5.022
beta_H[6,3] 8.392 1.181 6.668 8.120 10.878
beta_H[7,3] -2.558 0.623 -3.803 -2.548 -1.392
beta_H[8,3] 5.307 0.504 4.705 5.229 6.351
beta_H[9,3] -2.730 0.713 -4.106 -2.743 -1.340
beta_H[10,3] 8.545 0.270 8.013 8.544 9.082
beta_H[11,3] 8.519 0.287 7.904 8.539 9.039
beta_H[12,3] 5.269 0.334 4.491 5.313 5.779
beta_H[13,3] 8.859 0.174 8.501 8.862 9.191
beta_H[14,3] 5.790 0.273 5.187 5.814 6.269
beta_H[15,3] 10.276 0.294 9.708 10.270 10.858
beta_H[16,3] 6.030 0.618 4.742 6.060 7.116
beta_H[1,4] 8.102 0.188 7.711 8.110 8.450
beta_H[2,4] 9.850 0.129 9.584 9.857 10.079
beta_H[3,4] 9.952 0.153 9.626 9.963 10.217
beta_H[4,4] 11.819 0.464 10.865 11.835 12.707
beta_H[5,4] 5.555 0.724 4.349 5.485 7.223
beta_H[6,4] 7.133 0.914 5.247 7.324 8.468
beta_H[7,4] 8.542 0.348 7.821 8.542 9.207
beta_H[8,4] 6.727 0.262 6.219 6.734 7.199
beta_H[9,4] 7.438 0.454 6.518 7.445 8.327
beta_H[10,4] 7.525 0.216 7.117 7.518 7.959
beta_H[11,4] 9.515 0.187 9.147 9.513 9.891
beta_H[12,4] 7.260 0.196 6.886 7.256 7.662
beta_H[13,4] 9.071 0.138 8.797 9.076 9.324
beta_H[14,4] 7.843 0.214 7.443 7.835 8.284
beta_H[15,4] 9.477 0.220 9.050 9.477 9.893
beta_H[16,4] 9.399 0.242 8.948 9.394 9.883
beta_H[1,5] 8.907 0.152 8.590 8.907 9.202
beta_H[2,5] 10.676 0.100 10.483 10.675 10.877
beta_H[3,5] 10.755 0.161 10.460 10.748 11.088
beta_H[4,5] 8.344 0.487 7.406 8.326 9.339
beta_H[5,5] 5.194 0.598 3.672 5.259 6.156
beta_H[6,5] 8.963 0.628 8.014 8.856 10.324
beta_H[7,5] 6.587 0.347 5.926 6.582 7.306
beta_H[8,5] 8.249 0.201 7.919 8.234 8.641
beta_H[9,5] 8.197 0.468 7.251 8.201 9.116
beta_H[10,5] 10.039 0.222 9.604 10.039 10.476
beta_H[11,5] 11.584 0.218 11.159 11.583 12.004
beta_H[12,5] 8.611 0.186 8.223 8.613 8.974
beta_H[13,5] 10.142 0.130 9.886 10.139 10.401
beta_H[14,5] 9.235 0.213 8.831 9.230 9.675
beta_H[15,5] 11.302 0.236 10.833 11.303 11.761
beta_H[16,5] 9.844 0.182 9.463 9.853 10.177
beta_H[1,6] 10.114 0.231 9.736 10.084 10.665
beta_H[2,6] 11.473 0.118 11.244 11.473 11.702
beta_H[3,6] 10.618 0.156 10.287 10.626 10.903
beta_H[4,6] 12.908 0.883 11.142 12.928 14.558
beta_H[5,6] 5.654 0.607 4.592 5.624 6.936
beta_H[6,6] 8.715 0.578 7.423 8.749 9.788
beta_H[7,6] 10.537 0.590 9.341 10.552 11.705
beta_H[8,6] 9.527 0.253 9.048 9.538 9.954
beta_H[9,6] 8.278 0.785 6.738 8.268 9.847
beta_H[10,6] 9.152 0.346 8.418 9.178 9.764
beta_H[11,6] 10.711 0.360 9.947 10.736 11.340
beta_H[12,6] 9.377 0.237 8.939 9.367 9.869
beta_H[13,6] 11.092 0.148 10.805 11.090 11.389
beta_H[14,6] 9.918 0.279 9.338 9.923 10.458
beta_H[15,6] 10.571 0.422 9.701 10.573 11.408
beta_H[16,6] 10.605 0.248 10.100 10.609 11.089
beta_H[1,7] 10.206 1.204 7.164 10.500 11.877
beta_H[2,7] 12.104 0.532 10.919 12.141 13.070
beta_H[3,7] 10.295 0.665 8.854 10.351 11.403
beta_H[4,7] 1.468 4.636 -7.299 1.352 10.630
beta_H[5,7] 5.865 1.908 2.174 5.912 9.702
beta_H[6,7] 7.023 3.213 -1.000 8.116 10.896
beta_H[7,7] 7.694 2.948 2.037 7.668 13.544
beta_H[8,7] 10.863 0.968 9.196 10.859 12.481
beta_H[9,7] 5.070 3.991 -3.007 5.129 12.720
beta_H[10,7] 9.994 1.613 7.228 9.833 13.550
beta_H[11,7] 11.338 1.793 8.160 11.165 15.334
beta_H[12,7] 9.762 0.999 7.596 9.849 11.396
beta_H[13,7] 11.720 0.708 10.052 11.789 12.924
beta_H[14,7] 10.633 0.892 8.799 10.678 12.250
beta_H[15,7] 13.862 2.204 9.554 13.864 18.267
beta_H[16,7] 12.738 1.362 10.353 12.610 15.718
beta0_H[1] 9.149 17.963 -26.037 8.939 45.422
beta0_H[2] 10.524 8.076 -4.773 10.545 27.092
beta0_H[3] 9.631 10.734 -10.702 9.734 30.369
beta0_H[4] 9.475 176.716 -345.977 7.494 364.342
beta0_H[5] 4.721 27.708 -46.460 4.318 57.761
beta0_H[6] 8.192 52.992 -98.212 7.771 129.531
beta0_H[7] 5.820 143.553 -288.824 6.375 308.591
beta0_H[8] 6.196 24.583 -19.467 6.430 31.705
beta0_H[9] 5.741 123.768 -249.177 5.480 251.950
beta0_H[10] 8.333 34.744 -63.328 8.964 76.485
beta0_H[11] 9.531 50.108 -100.927 9.882 114.943
beta0_H[12] 6.522 13.098 -18.760 6.337 32.405
beta0_H[13] 9.781 10.390 -11.381 9.762 29.457
beta0_H[14] 7.102 11.019 -15.804 7.037 27.830
beta0_H[15] 9.147 116.228 -225.346 9.307 252.444
beta0_H[16] 7.239 28.399 -51.772 7.413 64.850